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Abstract 

We report recent progress on theoretical investigations of quarkonia at finite tem- 
perature. We discuss medium modification of charmonia and bottomonia from a view- 
point of local operators and point out that while charmonia are sensitive to the decon- 
finement transition bottomonia will be modified at much higher temperatures. 
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§1. Introduction 

Properties of heavy quarkonia have been extensively studied since it was pointed out 
that they provide information on deconfinement transition of QCDP'EJ While expected 
suppressions of the quarkonia have been measured in heavy ion experiments,® interpretation 
of these data is not so straightforward, because of not only the complexity of the collision 
processes but also the fact that spectral properties of the heavy quarkonia are not well 
understood yet. In this report, we focus on recent development of theoretical understanding 
of the quarkonium states at finite temperature. 

§2. Theoretical development 

Lattice QCD provides a unique non-perturbative first principle approach. A direct study 
of the quarkonium spectral function p(u, T) is possible with the help of the maximum entropy 
method (MEM) which enables us to invert a current correlation function at imaginary time 
G(t, T) via a dispersion relation.® 

ri t\ I A cosher -l/(2r))] 

G(T ' T) = J dU sinh[c/(2T)] (2 ' 1} 

Even in recent calculations with bigger lattice sizes,® however, quantitative information 
seems hard to be extracted, presumably due to limited temporal lattice size at high tem- 
perature and a small number of data points. The existence of the spectral peak does not 
necessary mean survival of a quarkoniunP and substantial spectral modification does not 
contradict with the behavior of G(t,T)P A rather promising direction seems to determine 
an interquark potential V(r, T) containing both real and imaginary parts for a Schodinger 
equation by lattice QCDP Indeed, existence of the imaginary part in the potential was 
pointed out by Laine et al. within a resummed perturbation theoryP Recently significant 
progress has been made with an effective field theory framework for heavy quark bound 
states.^* Though this approach assumes a hierarchy of the energy scale which becomes 
complicated at finite temperature due to newly introduced medium energy scales, analyti- 
cally results have been presented in some cases to give an insight of possible mechanisms of 
the in-medium modification of a quarkonium.^ These analyses indicate that a dominant 
in-medium effect on quarkonium at experimentally accessible temperatures could be a singlet 
to octet breakup process by gluons. On the other hand, an estimation of the energy density 
achieved in heavy ion collisions at RHIC energies and subsequent hydrodynamic evolution 
lead to a lifetime of the deconfined phase long enough to melt quarkonia if J 1 > 50 MeV.^ 
This means that we need a theoretical estimation of the width at those temperatures, likely 
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in the strongly interacting regime.^ 1 One of the promising approach for this purpose is 
relating local operators to a medium modification of a quarkonium via an operator product 
expansion (OPE). We will give a basic concept and a recent result below. 

§3. Local operator approach for quarkonia 

The interaction between a heavy quarkonium and soft gluons was first formulated by 
Peskin.^ The key concept is a separation scale, which is binding energy e in the case 
of a heavy quarkonium. Regarding an exchange momentum k larger than e as hard scale 
while the other soft one, we may express a matrix element via OPE as 

where C; and (Oi) stand for the Wilson coefficients responsible for the hard scale k > e and 
expectation value of local operators for the soft scale, respectively (see Fig. [1]). The leading 
order contribution is given by dimension four gluon condensate. At lower temperature than 
the seperation scale, one may impose all the medium effect on the change of the expectation 
value of the operators, which can be extracted from lattice calculations. Taking the real 
part of the matrix element immediately leads to a formula of a mass shift, which is the 
second order Stark effect in QCD, as Aijiq-q = —^--f (^AE 2 ^ where a® is the Bohr 
radius of the Coulombic bound state and (^AE 2 ^ is the temperature dependent part of 
the chromoelectric condensate. The expectation value of the dimension four gluon operators 
can be extracted from pressure p, energy density e and effective coupling constant a^ ff (T) 
obtained in lattice gauge theory as 

= TT^N- M °W + !^W) (3-1) 

where M and M 2 correspond to the gluonic part of the QCD trace anomaly, e — 3p and 
that of enthalpy density, e + p, respectively and shown in Fig. [TJ The rapid change of the 
energy density results in an abrupt increase of i^E 2 } thus downward mass shift of a heavy 
quarkonium in the vicnity of the phase transition. 

The above framework with separation scale e is applicable for a deeply bounded heavy 
quark-antiquark system. In reality, this condition might be questionable for charmonium at 
high temperature. One can turn to the current correlation function II (q 2 ) = J d i xe iq x (Tj(x)j(O)) 
in terms of OPE by going to deep Euclidean region in momentum space, q 2 = —Q 2 . Large 
negative q 2 enables us to compute II(q 2 ) in a perturbative manner with a pointlike current 
such as j M (x) = c(x)7 M c(x) and makes convergence property better than the above case. The 
relation to the physical quarkonia, q 2 = m%n, is kept through the dispersion relation. After 
the Borel transformation, which optimizes the dispersion relation such that the integral over 
the energy is dominated by the lowest resonance, the dispersion relation for the transformed 
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T/T c 

Fig. 1. Left: Schematic diagram for OPE. Right: Temperature dependence of Mq/T 4 and M2/T 
extracted from pure SU(3) lattice gauge theory.® 

correlator M(y) reads 

M{u) = J dx 2 e- ux2 p(2m Q x,T) (3-2) 

where v = Atuq/M 2 and M is the so-called Borel mass parameter. QCD sum rules® gives a 
systematic framework to extract spectral properties from the current correlation function and 
dispersion relation also at finite temperature by introducing medium dependent condensates 
unless the typical energy scale of the medium exceeds the separation scaled With a Breit- 
Wigner type pole ansatz for the physical spectral function, one can derive a constraint on 
the spectral change of a quarkonium at finite temperature.®'^ 1 Recently, Gubler and Oka 
proposed to apply MEM to QCD sum rules.® In this method, we do not have to assume a 
specfic functional form on the spectral function. Futhermore, compared to lattice calculations 
based on Eq. (12-ip . we can take as many points as possible and the dispersion relation does 
not have a temperature dependence other than the spectral function itself. Figure [2] shows 
spectral functions of J/ip (left) 211 and T (right)® obtained by the QCD sum rule+MEM 
approach. Although resolution of the width in the lowest peak is not sufficient, one sees 
how the peak dissolves as temperature increases. The drastic change around T c seen in the 
case of J ftp is consistent with previous sum rule calculations,®'^ while T hardly reflects 
the phase transition but exhibits sizable modification above 2T C . 

§4. Summary and outlook 

Several approaches have revealed that a heavy quarkonium has a width at finite temper- 
ature induced by the medium effect, especially due to gluonic dissociations. Utilizing OPE, 
we have shown that charmonia are sensitive to the phase transition while bottomonium 
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Fig. 2. Spectral function obtained from QCD sum rules with maximum entropy method. Left: 
J/ip, Right: T. 



are modified at much higher temperatures. Since bottomonia will be more appropriate for 
treatment with effective field theories, we expect theoretical attempts could meet together 
to unveil the interaction of the heavy quarkonium in the deconfined medium, which can be 
explored in heavy ion colllisions at the LHC energies. 
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